Bursting: when a cusp and a pitchfork interact 
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We present an experimental and theoretical study of an unusual bursting mechanism in a two- 
mode semiconductor laser with single-mode optical injection. By tuning the strength and frequency 
of the injected light we find a transition from purely single-mode intensity oscillations to bursting 
in the intensity of the uninjected mode. We explain this phenomenon on the basis of a simple 
two-dimensional dynamical system, and show that the bursting in our experiment is organised by a 
cusp-pitchfork bifurcation of limit cycles. 
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Bursting is a striking feature in a wide range of physi- 
cal and biological systems. It is characterised by periods 
of quiescence that are interspersed with periods of activ- 
ity. Well known systems that exhibit bursting include 
the Hodgkin-Huxley model [l] |2j and Taylor-Couette 
flows [21 H] . Understanding the mechanisms responsible 
for bursting is a fundamental problem for dynamical sys- 
tems theory and the physical sciences in general [S^i^. 

Systems with invariant manifolds are abundant in 
practical applications. Interesting dynamical phenomena 
are associated with the loss of transverse stability of an 
attractor contained in an invariant manifold [8, 9j. For 
example, a blowout bifurcation leads to a form of burst- 
ing known as on-off intermit tency. Here, the quiescent 
phase corresponds to chaotic dynamics conflned to the 
invariant manifold. We recently verifled experimentally 
that on-off intermittency occurs in an optically injected 
two-mode laser system [10]. 

In this letter we study theoretically and conflrm ex- 
perimentally a new mechanism that leads to bursting in 
the optically injected laser system. The mechanism is a 
global heteroclinic bifurcation which is organised by an 
interaction between a cusp and a pitchfork bifurcation of 
limit cycles. In contrast to on-off intermittency, the qui- 
escent phase is governed by a nearby saddle-node bifur- 
cation within the invariant manifold, and a transversally 
unstable limit cycle. The essential dynamics responsible 
for this bursting are captured by simple two-dimensional 
normal form equations for the cusp-pitchfork bifurcation 
of equilibria [TT] . 

Our experiment employs a uni-directional coupling 
scheme as shown in Fig. [l] The master laser is a narrow 
linewidth (100 kHz) tunable laser. The two- mode semi- 
conductor laser is operated at twice threshold, where, in 
the absence of optical injection, weak coupling between 
the two modes leads to stable two-mode lasing (see the 
red optical spectra in Fig. [T]). Details of the design and 
lasing characteristics of this device are given in Ref. [12 . 
Light emitted from the master laser is injected into the 
long wavelength mode of the two-mode laser. Depend- 
ing on the strength and frequency of the injected light. 
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FIG. 1: (Colour online) (a) Single- mode light from the tun- 
able master laser is optically injected into the two-mode semi- 
conductor laser, (b) Optical spectrum of the free running two- 
mode laser. The light is injected close to the long wavelength 
mode as indicated by the arrow. 
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FIG. 2: Bursting time series of the uninjected mode obtained 
(a)-(c) experimentally and (d)-(f) theoretically. Acj and K 
take fixed values indicated by the white horizontal line in 
Fig. [3] From top to bottom the value of K is decreased. 



drastic changes in the output of the two-mode laser can 
occur. These changes are visible in the time series of 
the individual modes which we recorded using high-speed 
photodiodes. 

The bursting we are interested in occurs in a transi- 
tion between single-mode dynamics, where the uninjected 
mode is suppressed and the intensity of the injected mode 
oscillates periodically, and two-mode dynamics, where 
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both modes emit light with significant intensity. Figure |2] 
shows time series of the uninjected mode during this tran- 
sition for parameter values indicated by the white hori- 
zontal line in Fig.jS] A dramatic increase in the frequency 
of bursts is seen as the injection strength is decreased. 

A simple four-dimensional set of ODEs reproduces 
this bursting behaviour (Fig. [2|d)-(f)), along with other 
experimentally observed phenomena [T3| for the setup 
shown in Fig. [l] The normalised, time-dependent vari- 
ables in this model are: the magnitude of the uninjected 
mode, \Ei{t)\^ the complex field of the injected mode, 
£^2(^)7 and the excess carrier density, N{t). In dimen- 
sionless form it is given by \TS\ : 



d\Ei\ 
dt 
dE2 
dt 

rr^dN 

T-r- 
dt 



= -((2iV + l)5i-l)|i?i|, 



{{2N +I)g2-l){l + ia) -iAuj 



E2 + K, 



P-N-{l + 2N)[g,\Eif+g2\E2f], (1) 



where time is measured in units of the photon lifetime, 
^ph = 1/(9.8 X 10^^) s. The bifurcation parameters are 
the normalised injection strength, and the normalised 
frequency difference (or detuning) between the injected 
mode and the injected light, Auj. Further parameters are 
the coupling, a, between the magnitude and phase of E2; 
the product of carrier lifetime and the cavity decay rate, 
T; and the normalised pump current, P. The two modes 
are coupled through N{t) via the nonlinear modal gain 

g^=[l + e{\E,\'+/3\E2\^)]-' , 
g2=[l + e{/3\E,\^ + \E2f)]-' , 

where e is the self saturation of a mode, and e/3 is the 
cross saturation between modes. To agree with the ex- 
periments in Ref. [13] we fix e = 0.01, /3 = 2/3, a = 2.6, 
T = 800, and P = 0.5 (twice threshold). 

The phase space of system ([T]) features a three- 
dimensional sub-manifold, 



M = { {\Ei\,Re{E2),lm{E2),N) e 



\Ei\=0 }, 



which is invariant under the time evolution of ([T]). We 
refer to A4 as the single-mode manifold since the dynam- 
ics on A4 reduce to that of the well-studied single-mode 
injection problem [Ml [15]. System ([T]) has a formal Z2 
symmetry due to the transformation |£^i| — Only 
the positive magnitudes are relevant but the Z2 symme- 
try alters the generic set of bifurcations for system ([T]) . 

Figure [3] contains two-parameter bifurcations diagrams 
in the (K, Aa;)-plane. In the experimental bifurcation 
diagram (Fig. [sj^a)) regions of the parameter plane are 
coloured to represent the number of dominant peaks 
in the power spectrum of the uninjected mode for fre- 
quencies less than 1.5 times the relaxation oscillation. 
Figure [sj^b) contains a partial bifurcation diagram for 
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FIG. 3: (Colour online) Two parameter bifurcation diagrams. 

(a) experimental: colours represent the number of dominant 
peaks in the power spectrum of the uninjected mode, dark- 
grey = 0, light-grey = 1, green= 2, and yellow-red= 3-9. 

(b) numerics: bifurcations of limit cycles are saddle-node 
(SL), pitchfork (PL), period doubling (PD), torus (T), and 
saddle-node-pitchfork (SPL). Dashed lines indicate bifurca- 
tions solely of single-mode objects. Background colours in- 
dicate regions of periodic (grey), quasi-periodic (blue), and 
chaotic oscillations (yellow- red ) . Dark shading indicates that 
|^i(t)|=0. 



system ([T]). Curves are bifurcations continued using 
AUTO [16 and the coloured background represents dif- 
ferent attractor types characterised by their Lyapunov 
spectrum [17 . Additionally, parameter regions are 
shaded if the corresponding attractor is contained in the 
single-mode manifold M. (i.e. \Ei{t)\ = 0). The qualita- 
tive agreement between experiment and theory is clear. 
The dark-grey region in Fig. |3ja) corresponds to the un- 
injected mode being off. Leaving the dark- grey region, 
e.g. by increasing Acj, results in the uninjected mode 
turning on. Exactly how the uninjected mode turns 
on depends on where the dark-grey region's boundary 
is crossed. On the boundary between dark- and light- 
grey are supercritical pitchfork of limit cycle bifurcations 
(Fig.jsjb)). Crossing this boundary causes the uninjected 
mode to gradually turn on. There are two boundaries 
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where the uninjected mode turns on in a non-trivial man- 
ner. Both of these boundaries are associated with burst- 
ing dynamics. In the experimental bifurcation diagram 
(Fig. [3]^ a)), these are the boundaries between dark-grey 
and yellow/red regions. On-off inter mittency was found 
along the boundary at low Auj [10 . In this letter we focus 
on the second boundary at larger Auj. On this boundary 
there are two saddle-node pitchfork of limit cycle (SPL) 
bifurcations. These are formed by tangencies between 
a curve of saddle-node of limit cycle bifurcations and a 
curve of pitchfork of limit cycle bifurcations (Fig. |3jb)). 
Connecting the two SPL points is a curve of supercritical 
torus bifurcations. As we will show, this arrangement of 
bifurcations, which we refer to as the torus bubble, is due 
to a nearby interaction between a cusp and a pitchfork 
bifurcation of limit cycles. 

In order to gain analytic insight into the bifurcation 
structure, we first study a simplified system where the 
limit cycles are replaced by equilibria. We demonstrate 
that a large part of the phenomena observed in Fig. 3, 
including the bursting shown in Fig. 2, are a consequence 
of the unfolding of a cusp-pitchfork bifurcation [11 . Sim- 
ilar to Ref. [11 we choose a simple two-dimensional nor- 
mal form for a system which is close to a cusp-pitchfork 
bifurcation, 

r = r(^i + z) , 

i = -6 - ^3Z -r^-z\ (2) 

The bifurcation parameters are ^1, ^2, ^3 and we note 
that the phase space (r, z) of system ([2| has an invariant 
manifold defined by r = 0, analogous to the single-mode 
manifold M. of system ([T]). The local bifurcation struc- 
ture of system ([2| is shown in Fig. [4j Saddle-node bifur- 
cations take place on the (blue) surface that resembles an 
open book. For parameter values on the outside of the 
open book (^3 > —3(^2/2)^^^), system ^ has one equi- 
librium in the invariant manifold. For parameters values 
inside the open book (^3 < -3(^2/2)^^^), system ^ has 
three equilibria in the invariant manifold. The (purple) 
folded surface in Fig. [4] corresponds to pitchfork bifur- 
cations. Below the pitchfork surface system ([2| has no 
equilibria off the invariant manifold. Crossing through 
the pitchfork surface causes two symmetrically related 
equilibria to bifurcate out of the invariant manifold; one 
with r > and the other with r < 0. The remaining 
(brown) surface corresponds to Hopf bifurcations where 
each of the symmetrically related equilibria bifurcate to 
a limit cycle. The three surfaces meet at the black curve 
which comprises two branches of saddle-node-pitchfork 
bifurcations that emanate from the cusp-pitchfork bifur- 
cation at the origin. 

The transparent surface in Fig.[4]indicates a slice of pa- 
rameter space. Intersections of the bifurcations with this 
surface are projected down to the bottom of the figure. 
From these projections we see the same arrangement of 




FIG. 4: (Colour online) Unfolding of a cusp-pitchfork bifur- 
cation. Saddle- node surface (blue), pitchfork surface (pur- 
ple), Hopf surface (brown), and saddle-node-pitchfork curve 
(black). The bifurcation curves on the bottom are projections 
of the intersection of the transparent surface with the bifurca- 
tions. Bifurcations of equilibria are saddle- node (S), pitchfork 
(P), Hopf (H), and saddle-node pitchfork (SP). 



bifurcations as in the torus bubble (Fig. [3|, but for equi- 
libria instead of limit cycles. 

Figure [5] contains a sketch of the phase portraits 
around the saddle-node pitchfork (SP) point closest to 
the cusp-pitchfork bifurcation. We focus on the curve 
of global saddle-node heteroclinic (SHET) bifurcations. 
Along the SHET bifurcation curve, system ([2| has a 
transversally stable, nonhyperbolic equilibrium; and a 
heteroclinic cycle connecting the nonhyperbolic equilib- 
rium to a transversally unstable saddle equilibrium. For 
parameter values on one side of the SHET curve (re- 
gion (ii)) the nonhyperbolic equilibrium splits into two 
equilibria and the heteroclinic cycle is destroyed. For 
parameter values on the other side of the SHET curve 
(region (iii)) the nonhyperbolic equilibrium is destroyed 
and the heteroclinic cycle bifurcates to a stable limit cy- 
cle. In region (iii), close to the SHET curve, the ghost 
of the saddle-node bifurcation results in a slow region of 
phase space. While a trajectory of system ([2| is in this 
slow region it is attracted towards the invariant manifold 
(quiescent phase). Once through the slow region the tra- 
jectory shadows the dynamics in the invariant manifold, 
and is thus attracted toward a saddle equilibrium which 
is stable within the invariant manifold but transversally 
unstable. This equilibrium repels the trajectory away 
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FIG. 5: (Colour online) Bifurcation diagram for the saddle- 
node-pitchfork bifurcation, (i)-(v) phase portraits in the 
(r, z) plane. The vertical lines corresponds to the in- 
variant manifold (r = 0). Closed triangle=stable equi- 
librium, open triangle=saddle-equilibrium, and the closed 
square=nonhyperbolic equilibrium. 
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FIG. 6: (a)-(c) periodic bursting for K = 0.0092473 and 
Acj = 3.088 GHz. (d)-(f) chaotic bursting for K = 0.0085 
and Auj = 3.5126 GHz. 



from the invariant manifold (active phase) and the tra- 
jectory is then globally reinjected via the limit cycle to 
the slow region of phase space for the process to start 
over again. 

The mechanism responsible for bursting in the two- 
mode laser (Fig. |2| is analogous to that described in the 
previous paragraph but for limit cycles instead of equi- 
libria. We illustrate this in Fig. [6]^a)-(c). Figure [gJ a) 
shows bursting in the time series of the uninjected mode. 
The same plot on a log scale (Fig.|6|b)) reveals an almost 
linear trend towards the single- mode manifold A4. The 
slope of this linear trend is related to the negative trans- 
verse Lyapunov exponent of the nearby saddle-node limit 
cycle. Similarly, there is a phase of linear trend away 
from A4. The slope of this trend is related to the pos- 
itive transverse Lyapunov exponent of the saddle-cycle 



involved in the nearby SHET bifurcation. Moreover, two 
distinct periodic-like oscillations can be seen in the time 
series of the laser's population inversion (Fig. [gJc)). The 
time-scale of the quiescent phase is strongly influenced 
by the slow region of phase space. Since the slow region 
is caused by the ghost of a saddle-node of limit cycle 
we can estimate that the quiescent phase scales as /i~ 2 ^ 
where /i is the distance from the saddle-node of limit cy- 
cle bifurcation. We were able to experimentally verify 
this scaling. 

While the bifurcations responsible for the bursting in 
Fig. [6]^a)-(c) are similar to those in system the cusp- 
pitchfork bifurcation in the two-mode laser (Fig. [3| is of 
limit cycles, not equilibria. The extra rotation can com- 
plicate the bifurcation structure. For example, the stable 
torus could break up before the saddle- node of limit cycle 
bifurcation, and any curves of global bifurcations could 
be replaced by cones of heteroclinic tangles [18l|T9]. Fur- 
thermore, additional bifurcations can occur away from 
the cusp-pitchfork bifurcation. These complications can 
lead to different types of bursting. Such an example is 
shown in Fig. [6]^d)-(f). Here, the transversally unstable 
saddle-cycle is replaced by a chaotic saddle whose un- 
stable direction is transverse to M. The bursts are not 
(quasi-)periodic and the amplitude of the bursts varies 
from burst to burst (Fig. [oj^d)). The same plot on a log 
scale (Fig. [6]^e)) reveals phases of a linear trend towards 
A1, again this is associated with the negative transverse 
Lyapunov exponent of the nearby saddle-node limit cy- 
cle. The diffusive nature of the growth away from M. 
is strongly indicative of a transversally unstable (to A4) 
chaotic saddle. When a trajectory is very close to M, a 
longer signature of the chaotic saddle is observed before 
bursting away from M (Fig. |6|f)). 

In summary, we identified the cusp-pitchfork bifurca- 
tion of limit cycles as a mechanism for experimentally ob- 
served bursting in a two-mode semiconductor laser with 
optical injection. The backbone of the bifurcation struc- 
ture responsible for this bursting was explained by sim- 
ple two-dimensional normal form equations for the cusp- 
pitchfork bifurcation of equilibria. 
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